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Abstract 

We present a short review of the status of the dark matter problem. In 
particular we show that one of the best motivated candidate for the dark 
matter is the neutralino, a supersymmetric particle. Finally we study 
the possibility to detect 7-rays coming from neutralino pair annihilations 
that take place in our Galactic Center with the upcoming satellite detector 
GLAST. 

1 Introduction 

The dark matter problem is one of the most fascinating and intriguing issue 
in the astroparticle physics. This subject requires theoretical concepts and 
experimental techniques coming from cosmology, astrophysics and fundamental 
particle physics. 

By definition dark matter is a kind of non luminous matter. Its existence is 
required by a bunch of striking experimental evidences as well as by some com- 
pelling theoretical motivations. One of most interesting experimental evidences 
is given by the observations of the spiral galaxy rotation curves 3 . They can be 
determined well outside the luminous core of the galaxy, measuring the circular 
velocities of the neutral hydrogen clouds using the A = 21 cm lines emissions. 
The result when expressed as a function of the galactocentric distance R is that 
vc ~ const, for large R (in particular for R > Rium — 5Kpc) rather than the 
expected vc oc R^^/^. If we suppose that the first Newton's law still holds at 
the galactic scale we must conclude that M(< R) (x R where M(< R) is the 
total mass inside a shell of radius R. Hence there must be a "dark" contribution 
to the galactic mass M. 

Recent measurements ^ of the anisotropy of the Cosmic Microwave Back- 
ground (CMB) strongly indicates that the total matter contribution is flh^ ^ 0.3 
where h ^ 0.71 is the Hubble constant in units of 100 km s^^ Mpc^^. Hence 



1 



fim/i^ is much larger than the baryonic term njjh'^ 0.02. Necessarily the dom- 
inant matter component must be of a non baryonic form. It turns out that the 
best motivated candidates are some kind of Weak Interacting Massive Particles 
(WIMPs). Usually these type of candidates are stable particles that appear in 
various extension of the Standard Model (SM) of particle physics. 



2 Supersymmetric dark matter 

From the theoretical point of view the most interesting extension of the SM 
involves supersymmetry. Here we will focus in particular on iV = 1 supersym- 
metric extension of the SM with soft supersymmetry breaking terms. This is 
the so called Minimal Supersymmetric Standard Model (MSSM) P]. Usually 
the Lightest Supersymmetric Particle (LSP) that appear in the MSSM mass 
spectrum is a good CDM candidate It turns out that the LSP is often the 
lightest neutralino. The four neutralino gauge eigenstates are linear combination 
of the gaugino and higgsino fields: 

X^ = NaB + Ar,2VF° -I- NaH^a + N.^H^ (1) 

where i = 1,...4. The corresponding mass eigenstates can be obtained by 
diagonalizing the mass matrix: 



f Ml -CfjSwmz Sfsswmz \ 

M2 c/scwiTiz -SfjCwfnz 

-cpswmz C/scwmz -fi 

\ spswmz -spcwmz ^ / 



(2) 



that depends by various MSSM parameters, like the gaugino masses Mi and 
M2 (that are associated, in the lagrangian, to explicit supersymmetry breaking 
terms) or the fi parameter that appears in the Higgs sector [1]. The lowest 
eigenvalue of the mass matrix in equation |21 is exactly the LSP, until the so 
called R-parity is conserved The neutralino is a Majorana fermion, so it 
coincides with its own antiparticle, and it is a weak interacting particle. The 
R-parity conservation implies that the lightest neutralino cannot decay in SM 
particles and so it is a stable particle. For these reasons the neutralino is a good 
CDM candidate. Moreover a neutralino pair can self annihilate into various SM 
particles that we can detect. 



3 Supersymmetric dark matter with GLAST 

In the previous section we have seen that the neutralino is a good CDM candi- 
date. It is possible to study its properties through indirect detection of cosmic 
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rays coming from pair annihilations in the dark galactic halo. In particular here 
we focus on the cosmic 7-rays from the Galactic Center (GC). 

The EGRET data fr'om the GC |3 give a strong indication of an excess 
with respect of the "standard" model production of 7-rays fBl . In the standard 
picture, the diffuse 7-ray background is due to the following processes: 

• tt'^ production that promptly decays in 27 

p + X ^ .. ^ 7r° ^ 27, i/e + X ^ .. ^ 7r° 27 

where p and He are primary protons and heliums and X is the interstellar 
hydrogen or helium. 

• Bremsstrahlung 

• Inverse Compton 

We are only interested in the energy range above E > 1 GeV where the dominant 
component is by far the tt'^ production. In the subsequent analysis we consider 
only this component for the diffuse background. 

We want to explain the GC excess adding a neutralino induced component 
in the tt'^ production. The leading intermediate annihilation channels for the 
tt'^ production, through fragmentation and/or decay processes, are often bb, tt, 
W'^W~ and Z'^Z'^ (see figure^. This is valid not only for the neutralino but 
also for a generic Majorana fermion WIMP, as for such particles the s-wave 
annihilation rate into lighter fermions is suppressed by a factor m^/m^, where 
m/ is the final state fermion mass and is the WIMP mass. 

We can write the total diffuse 7-ray flux as a sum of two components: 

<^™^ = (/.;^^« + (/.^x = NkSb + ^x'^x (3) 
where the background flux cj)^^^ is splitted into two factors: 

[1 cm^sr) 

. , \ ,7 (4) 
(Icm-'sr-i) 4» 4!l"'{l = IS) 

where Em,{Ey) (measured in GeV~^s^^) is the local emissivity per hydrogen 
atom. The neutralino induced component N^(j>^ is given by 0: 

s;^ai\f ^ 

5/ 



ME,) = 3.74. 10^^" " , )('-^) y:^Bj 

V 10-26 cm^s-iy V J ^ dE ^ 



{JmAn^^J^^dn'Jiij') (5) 
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Figure 1: Relevant intermediate annihilation channels for a Majorana fermion 
WIMP (in particular for a neutralino). 



where Afl is the detector angular acceptance and where we have introduced the 
dimensionless function J that contains the dependence on the dark halo density 
profile: 
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f 



(6) 



A.5Kpc \0.3GeV/cm^y 
We have introduced the unknown normalization factor Nh in order to take into 
account our ignorance of the exact hydrogen column density, i.e. the interstellar 
medium, while the factor parametrize our ignorance of the dark matter halo 
model. Results from N-body simulations 0^13 give for for the dark matter 
density: 



p{r) - Po (^) 
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{ro/ar 



{p-l)/a 



1 + (r/a)" J 

where a, (3 and 7 are parameters that describe the halo profile (Isothermal 
sphere, NFW, Moore, ..), po is the WIMP density measured here and ro is the 
galactocentric distance. 

Having introduced the neutralino induced component the fit of the GC 
EGRET data greatly improves. In the context of a simplified toy model [7], 
valid for a generic Majorana fermion WIMP of mass and with a fixed dom- 
inant intermediate annihilation channel, the result is shown in figure |2] Our 
analysis shows that the EGRET data fit improves for small neutralino masses. 
In order to extract much more information about the GC EGRET excess we 
have to wait for the upcoming generation of telescope satellite, like GLAST [5]. 
With respect to EGRET, GLAST has a wider energy range, increased effective 
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Figure 2: Fit of the GC EGRET data assuming the the excess is due to a 
WIMP (of mass = 80 GeV) induced component with W^W~ dominant 
annihilation channel. The best fit parameters Nh and are indicated in the 
upper right corner. 

area and better energy and angular resolution. Hence we have studied what 
kind of data GLAST would collect from the GC, under the hypotesis that the 
excess, as mapped by EGRET, is due to WIMP annihilations. Relying on a 
simplified picture of the GLAST detector performances (energy resolution of 
about 10%, angular resolution of lO^'"^ sr and peak effective area of ~ 10^ cm^), 
we have computed the GC data set which will be obtained by GLAST in 2 years. 
We have superimposed the error bars associated to the statistical errors only 
for the chosen energy binning and with an angular acceptance of Ail = 10^"^ sr. 
The result, in this case, is shown in figure [31 

All the previous results are valid in the case of a generic WIMP. At the 
end we want to identify the neutralino as our WIMP candidate. To do so we 
must consider a "realistic" supersymmetric model, i.e. the MSSM. In order 
to reduce the number of free parameters of the MSSM it is possible different 
supersymmetry breaking scenarios in which the supersymmetry breaking terms 
derive from an underlying high energy theory. One of the most widely studied 
frameworks is the so-called minimal supergravity (mSUGRA) or constrained 
MSSM (cMSSM) [ni- In the mSUGRA model there is an important assumption 
about the universality of all the coupling constants at the grand unification scale 
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Figure 3: Data set collected by GLAST in 2 years in the case the EGRET GC 
excess is due to a WIMP induced flux. 



(GUT). With these constraints the number of free parameters is only flve: 

'^1/25 'mo, sign{^), Aq, tan(3 

where mo is the common scalar mass, mi/2 is the common gaugino mass and 
^0 is the proportionality factor between the supersymmetry breaking trilinear 
couplings and the Yukawa couplings, tan/3 denotes the ratio of the vacuum 
expectation values of the two neutral components of the SU(2) Higgs doublet, 
while the Higgs mixing /x is determined (up to a sign) by imposing the Electro- 
Weak Symmetry Breaking (EWSB) conditions at the weak scale. In this context 
the MSSM can be regarded as an effective low energy theory. The parameters at 
the weak energy scale are determined by the evolution of those at the unification 
scale, according to the renormalization group equations (RGEs). Fixing tan/3, 
Aq and sign{iJ,) wo have performed an accurate scan in the (toq, mi/2) mSUGRA 
parameter space, looking for the minimum normalization factor needed to be 
able to single out the neutralino annihilation signal with GLAST. In this case we 
have taken advantage of the better GLAST angular resolution (Af2 = 10~^ sr) 
with respect to EGRET. In this deeper analysis we have also computed the 
neutralino relic density: 

n^h^ = Uhlh^ (8) 

pc 

solving the Boltzmann equation that describes the time evolution of the neu- 
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Figure 4: Contour plots in the mSUGRA (mo, 7711/2) plane, for the value of 
the normalization factor N^, that allows the detection of the neutralino 7 ray 
signal, with GLAST. The light shaded region corresponds to the cosmologically 
favoured region where 0.1 < O^ft,^ < 1, while the dark shaded one corresponds 
to models that are excluded either by incorrect EWSB, LEP bounds violations 
or because the neutralino is not the LSP. 



tralino^ number density n^{t): 



This equation can be easily solved numerically including resonances, threshold 
effects and all possible coannihilation processes |13j. 

The results, for two possible choices of the parameters tan f3, Aq and sign{fi), 
are shown in figure 01 The contour plots, in the mSUGRA (mo,TOi/2) plane, 
show the normalization factor needed for GLAST detection, at 3a level, of 
the neutralino induced 7 ray signal. 

4 Conclusions 

We have performed an analysis of the possibility of an indirect detection of 
supersymmetric dark matter with cosmic 7-rays. The most important indication 
is that there is indeed room for supersymmetric dark matter in the already 
available 7-ray data from EGRET. Moreover the computation of the expected 7- 
ray flux coming from dark matter annihilation with upcoming detector GLAST 
shows that it will be possible to study in more details a supersymmetric signal. 

^ valid also for a generic WIMP 
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